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ABSTRACT 

If  a  magnetic  bit  in  a  thin  film  memory  array  is  driven  in  the  hard  direction,  the  changing  flux 
produces  eddy  currents  in  the  surrounding  conducting  surfaces.  It  has  been  pointed  out  by  Pahm 
and  others  that  the  eddy  field  produced  by  these  currents  olways  tends  ta  restore  the  bit  mognet- 
izatian  ta  its  original  direction.  If  this  eddy  field  were  strong  enough  at  the  end  of  the  ward 
drive  pulse,  the  effect  could  be  made  the  basis  af  a  scheme  far  nondestructive  readout  (NDRO). 

Same  experiments  were  performed  an  structures  where  the  eddy  field  was  strong  enough  ta  restore 
the  bit.  The  thin  film  plane  used  contained  50  X  50-mil  X  500-A-thicl<  Permal  lay  films  with  uni¬ 
axial  anisotropy.  The  word  line  was  completely  wrapped  around  the  plane.  Ward-1  ine-ta-ward- 
line  spacings  of  1 1 . 25  and  7. 25  mils,  and  ward-1  ine  widths  af  50  ,  75,  and  150  mils  were  used. 
The  eddy  field  wos  measured  at  50  ond  100  nsec  ofter  the  switching  aperotian  ond,  assuming  an 
exponential  eddy  field  decay,  a  time  constant  was  colculoted. 

Two  theoretical  approximations  were  considered.  In  bath,  the  eddy  currents  were  assumed  ta 
flaw  only  in  the  ward  lines.  In  one,  the  thickness  of  the  lines  was  assumed  negligible.  In  the 
other  approximation,  the  thickness  af  the  striplines  was  considered,  but  was  assumed  ta  be  af 
infinite  width.  The  finite  width,  thin  stripline  approximation  seems  ta  be  in  better  agreement 
with  the  measured  values  af  the  time  constant  which  varies  between  76  and  200  nsec. 


Accepted  far  the  Air  Farce 
Franklin  C.  Hudson,  Deputy  Chief 
Air  Farce  Lincoln  Laboratory  Office 


*This  report  is  based  an  a  thesis  af  the  same  title  submitted  ta  the  Department  af  Electrical 
Engineering  at  the  Massachusetts  Institute  af  Technology  an  19  August  1963,  in  partial  fulfill¬ 
ment  of  the  requirements  far  the  degree  af  Master  of  Science. 


TABLE  OF  CONTENTS 


Abstract  iii 

I.  Introduction  1 

11.  Magnetic  Thin  Film  Memories  2 

III.  Theory  3 

IV.  Experimental  Apparatus  15 

V.  Experimental  Procedure  and  Results  19 

VI.  Suggestions  for  Further  Study  23 

Appendix  A  27 

Appendix  B  28 


V 


NONDESTRUCTIVE  READOUT  SCHEME 
FOR  THIN  FILM  MEMORIES 


I.  INTRODUCTION 

In  recent  years,  much  attention  has  been  given  to  magnetic  thin  film  structures,  particularly 
for  memory  devices.  Most  of  these  devices  use  thin  film  bits  of  Permalloy,  vacuum  evaporated 
in  a  magnetic  field,  to  give  the  resultant  film  uniaxial  anisotropy.  A  few  small  memories  have 
been  manufactured  commercially. 

The  basic  advantage  of  magnetic  thin  films  is  that  they  switch  much  faster  than  ferrite  de¬ 
vices.  Because  of  the  thinness  of  the  film,  the  bit  tends  to  exist  in  a  single  domain.  Under  the 
influence  of  short,  powerful,  driving  fields,  the  entire  magnetic  spin  system  of  the  bit  rotates 
as  a  unit,  held  together  by  the  strong  exchange  interactions  between  the  atomic  spins.  There¬ 
fore  the  bit  switches,  for  fields  above  a  certain  threshold,  by  coherent  rotation  (a  fast  process 
which  can  be  complete  in  nanoseconds)  rather  than  by  wall  nucleation  and  motion,  a  much  slower 
process. 

It  has  been  found  possible  to  obtain  nondestructive  readout  (NDRO)  from  thin  films  by  limit- 
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ing  the  drive  fields  to  small  values.  Unfortunately,  the  signals  from  the  memory  array  are  then 
quite  small.  Also,  the  information  stored  by  the  bit  is  susceptible  to  processes  which  gradually 
destroy  the  information.  Recently,  another  technique  (suggested  by  Pohm  and  others)  involves 
passive  loading  of  the  thin  film  bit.  The  currents  induced  in  the  passive  loading  structures  when 
the  film  is  switched  produce  fields  which  restore  the  bit  to  its  original  state  at  the  end  of  the 
read  operation.  Eddy- NDRO  is  much  less  susceptible  to  gradual  information  destruction  than 
drive- limited  NDRO.  Furthermore,  much  larger  drive  fields  may  be  used  and,  therefore,  much 
larger  signals  may  be  obtained  from  the  films. 

An  understanding  of  the  details  of  eddy- NDRO  requires  a  thorough  analysis  of  the  interaction 
between  the  film  and  the  eddy  fields  produced  by  the  film.  Pohm  surrounds  his  bit  with  a  shorted 
copper  turn,  thus  making  the  analysis  quite  straightforward.  It  would  be  much  more  practical, 
however,  if  the  passive  loading  could  be  produced  by  the  striplines  used  for  driving  the  bit.  In 
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this  case,  the  analysis  would  be  much  more  complex.  Although  Smay  did  some  rough  calcula¬ 
tions  concerning  this,  it  was  found  that  a  much  more  precise  and  rigorous  procedure  was  needed 
to  adequately  handle  the  problem. 

The  v^ork  presented  in  this  report  has  been  directed,  basically,  toward  an  understanding  of 
some  of  the  details  of  the  interaction  between  the  eddy  currents  and  the  film.  Most  of  the  effort 
was  devoted  to  a  consideration  of  the  time  constant  of  the  eddy  current  decay.  Some  equations 
for  the  calculation  of  the  time  constants  were  developed  and  compared  with  some  experimental 
results. 
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II.  MAGNETIC  THIN  FILM  MEMORIES 


I 

The  formation  of  Bloch  walls  is  energetically  unfavorable  in  very  thin  magnetic  films. 
Therefore,  the  magnetic  bit  can  usually  be  considered  as  a  single  domain.  It  is  important  to 
note  that  under  suitable  conditions  the  film  can  break  up  into  a  multidomain  structure. 

If  the  film  is  evaporated  with  a  magnetic  field  parallel  to  the  plane  of  the  substrate,  the  re¬ 
sultant  film  will  have  an  axis  of  anisotropy  in  the  direction  in  which  the  original  field  had  been 
pointed.  The  magnetization  of  the  film  normally  points  along  the  anisotropy  axis  in  either  of 

the  two  possible  orientations.  As  the  magnetization  vector  is  rotated  away  from  the  anisotropy 
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axis,  in  the  plane  of  the  film,  the  energy  of  the  film  goes  as  k  sin  A,  where  A  is  the  angle 
between  the  magnetization  and  the  anisotropy  axis,  and  k  is  the  anisotropy  constant.  This  con¬ 
tribution  to  the  energy  of  the  system  is  called  the  anisotropy  energy.  It  gives  rise  to  a  torque, 
called  the  anisotropy  torque,  which  tends  to  rotate  the  magnetization  vector  back  into  the  anisot¬ 
ropy  axis  orientation  through  the  smaller  angle.  This  torque  is  zero  at  the  unstable  equilibrium 
at  A  =  90°. 

A  common  memory  structure  is  shown  in  Fig.  1.  With  no  field  applied,  the  magnetization  of 
the  bit  can  point  in  either  direction  along  the  easy  axis,  the  name  normally  given  to  the  anisotropy 
axis.  The  word  lines  produce  a  field,  in  the  hard  direction,  which  rotates  the  magnetization  of 
the  bit  until  it  points  in  the  hard  direction.  The  resultant  change  in  flux  linking  the  sense  line 
produces  the  output  voltage,  which  is  positive  or  negative  depending  upon  the  original  direction 
of  the  bit  magnetization.  Information  can  be  written  into  the  bit  by  first  rotating  the  magnetiza¬ 
tion  into  the  hard  direction.  While  the  bit  pulse  is  on,  the  word  pulse  is  removed  and  the  bit 
magnetization  rotates  back  toward  the  anisotropy  axis,  pointing  in  the  direction  chosen  by  the 
bit  field. 

Real  systems  have  certain  imperfections  which  must  be  discussed.  First,  the  direction  of 
the  anisotropy  axis  differs  slightly  from  point  to  point  through  the  film.  In  the  film  on  which 
most  of  the  experimentation  was  done,  for  example,  90  percent  of  the  magnetization  is  within 


Fig.  1.  Thin  film  memory  bit,  simplified  view. 
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an  anisotropy  axis  spread  of  4.6°.  To  obtain  good  writing,  the  field  in  the  easy  direction  must 
be  strong  enough  to  rotate  most  of  the  magnetization  to  either  one  side  or  the  other  of  the  A  =  90° 
point  so  that,  upon  the  removal  of  the  word  drive,  most  of  the  magnetization  will  fall  back  to 
the  easy  axis  pointing  in  the  same  direction.  The  spread  in  anisotropy  axis  is  called  dispersion. 

Also,  it  is  often  impossible  to  align  the  word  lines  so  that  they  are  precisely  parallel  to  the 
easy  axis  of  the  bit,  since  the  direction  of  the  easy  axis  varies  for  the  various  bits  in  the  plane 
(this  is  known  as  skew).  Although  the  skew  angle  is  small,  the  word  drive  is  quite  large.  This 
means  that  the  component  of  the  word  drive  along  the  easy  axis  may  be  appreciable.  For  good 
writing,  the  field  in  the  easy  direction  must  be  large  enough  to  overcome  both  skew  and 
dispersion. 

In  the  normal,  destructive  readout  (DRO)  mode,  the  bit  breaks  up  into  a  multidomain  struc¬ 
ture  at  the  end  of  the  word  pulse  on  a  read  operation.  The  part  of  the  magnetization  pointing  to 
the  left  of  the  hard  axis  (A  =  90°)  falls  back  to  the  left,  and  the  part  of  the  magnetization  pointing 
to  the  right  of  the  hard  axis  falls  back  so  that  it  points  to  the  right  along  the  easy  axis.  The  in¬ 
formation  read  out  is  then  rewritten  into  the  bit.  Actually,  the  bit  is  never  a  truly  one- domain 
structure,  because  there  are  usually  small  walls  formed  around  impurities  and  at  the  edges  of 
the  bit. 

The  film  can  be  operated  in  a  current- limited  NDRO  mode  by  keeping  the  word  current  too 
small  to  ever  bring  the  magnetization  vector  close  to  the  hard  axis.  This  restricts  us  to  small 
sense  line  signals.  Another  limitation  arises  from  the  fact  that  writing  a  bit  produces  an  easy 
direction  field  in  all  the  bits  along  the  bit  line,  including  those  not  being  written.  This  field  may 
cause  the  small  domains  in  the  unseleeted  bits  to  grow  until,  after  many  such  operations,  the 
information  in  the  unseleeted  bit  is  destroyed.  Oeeasional  reading  and  rewriting  of  the  bits,  as 
we  do  in  the  DRO  mode,  helps  to  alleviate  this  problem. 

To  operate  in  the  eddy- NDRO  mode,  the  eddy  field  (produced  when  the  bit  is  driven  into  the 
hard  direction)  must  be  strong  enough  to  overcome  skew  and  dispersion  and  restore  the  bit  to  its 
original  condition  wlien  the  word  pulse  is  removed.  In  the  time  between  the  rise  and  fall  of  the 
word  pulse,  the  eddy  field  will  be  rapidly  decaying;  therefore,  the  word  pulse  must  be  quite 
short.  However,  its  amplitude  may  be  made  large  enough  to  obtain  the  full,  DRO  sense  line 
signal. 

It  is  conceivable  that  the  eddy  field  might  even  be  able  to  regenerate  a  partially  broken- up 
bit.  If,  for  example,  the  eddy  field  produced  from  a  bit,  of  which  10  percent  had  broken  up, 
was  strong  enough  to  restore  91  percent  of  the  magnetization,  the  signal  would  gradually  im¬ 
prove  with  successive  NDRO  operations. 

III.  THEORY 

An  analysis  of  eddy- NDRO  naturally  falls  into  two  parts.  First,  we  must  find  out  how  the 
eddy  field  produced  by  the  bit  behaves  and,  second,  we  must  find  out  how  this  field  affects  tlie 
bit. 

The  rotational  time  constant  of  the  isolated  film  is  on  the  order  of  one  or  two  nanoseconds, 
so  it  can  be  assumed  that  the  film  will  respond  instantaneously  to  the  applied  fields.  The  applied 
fields,  however,  are  not  uniform  and,  therefore,  different  portions  of  the  film  spin  system  ex¬ 
perience  different  torques.  The  exchange  forces  are  fairly  large;  therefore,  we  assume  that  the 
film  magnetization  rotates  as  a  unit  under  an  average  torque  produced  by  an  average  field.  If 
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the  magnetization  attempted  to  rotate  out  of  the  plane  of  the  film,  we  would  have  free  poles  on 

the  upper  and  lower  surfaees  of  the  film.  This  would  produee  a  very  large  field  outside  the  bit. 

This  field  would  have  a  large  energy.  Therefore,  the  magnetization  veetor  is  constrained  to 

rotate  in  the  plane  of  the  film.  For  our  square  bits,  the  angle  of  the  magnetization  to  the  easy 

axis  is,  to  a  good  approximation,  the  equilibrium  angle  between  the  torques  due  to  the  applied 

fields,  ineluding  the  eddy  field,  and  the  torque  due  to  the  anisotropy. 

Suppose  the  magnetization  per  unit  volume  is  u.  Then,  the  torque  upon  the  spin  system 

when  the  magnetization  has  rotated  an  angle  A  from  the  easy  axis  is  (sec  Fig.  Z)  uH  cos  A  — 

Z  ^ 

ull^  sin  A.  The  anisotropy  energy  is  k  sin  A,  where  k  is  the  anisotropy  eonstant.  Therefore 
the  torque  is  — k  sinZA;  M  =  u{l  eos  A  +  1  sin  A).  As  a  result,  the  value  of  dM  is  u(l  eos  A  — 

/  ^  y  y 

1^  sin  A)  (dA/dt)  dt.  The  eddy  field  at  time  t  due  to  a  ehsinge  in  the  magnetization  at  time  p 

■1^  sin  A)  (dA/dp)  dp. 


will  be  some  function  of  t  —  p;  dH 

i  '  pj 


-BdMf(t  -  p)  -  -Bf(t  -  p)  u(l^  eos  A 


The  total  torque  T  due  to  the  eddy  currents  at  time  t  is 


T 


—  Bu^  |eos  A  ^  (dA/dp)  eos  A  f(t  —  p)  dp  +  sin  A  ^  (dA/ dp)  sin  A  f(t  —  p)  dpj 


We  ean  integrate  by  parts  to  put  the  equation  in  a  more  usable  form.  For  example,  if  we  are 
considering  a  turn-on  problem,  where  A  at  time  zero  is  zero 


T  =-Bu 


eos  A  |sinA  — y  sin  A(df/dp)  dp 
r  eos  A  (df/ dp)  dp 


+  sinA  I— cos  A  +  f(t) 


(1) 


We  chose  B  so  that  f(0)  is  one.  If  the  drive  is  a  word  pulse  in  the  y- direction 

uH  cos  A  =  — T  +  k  sin  ZA  ;  (Z) 

y 

f(t)  is  probably  quite  close  in  form  to  that  of  a  decaying  exponential.  Once  we  find  it  (and  B),  it 
is  quite  simple  to  integrate  the  preceding  equations  to  find  the  dependence  of  A  on  time. 

If  the  word  drive  is  strong  enough  to  rotate  the  bit  magnetization  into  the  hard  direction  in 
a  time  which  is  much  less  than  the  decay  time  of  f{t),  we  need  not  do  the  integrals  at  all.  We 
merely  need  to  find  the  eddy  field  produced  by  reducing  the  x- component  of  M  suddenly  to  zero 
and  increasing  the  y- component  of  M  suddenly  to  u.  We  are  basically  interested  in  the  changes 
in  M  and  H  in  the  easy  direction. 


Fig.  2.  Bit  coordinate  axes. 


4 


It  is  now  necessary  to  find  the  eddy  field  produced  by  a  change  in  the  film  magnetization. 
First,  the  normal  modes  of  deeay  will  be  found.  A  normal  mode  is  a  field  distribution  whose 
time  dependence  is  completely  contained  in  a  multiplicative  factor  of  exp[— st]  .  We  know  that, 
immediately  after  the  magnetization  change,  an  eddy  current  will  flow  to  keep  the  field  in  the 
copper  constant.  This  transient  current,  and  the  transient  field  it  produces,  will  gradually  go 
to  zero  as  the  field  pattern  goes  to  its  new  steady  state.  We  ean  write  the  transient  field  as  the 
sum  of  normal  modes,  and  thereby  calculate  the  behavior  of  the  field  with  time. 

Under  the  assumption  of  exponential  deeay,  Maxwell's  equations  outside  the  film  become 

V  X  E  =  -  H 
c 

V  •  E  =  —  47rp 

c 

Wc  have  ignored  the  displacement  current  compared  with  the  conduction  current.  Where  p  is 
the  I'esistivity  of  the  copper,  E  =  pj.  Therefore,  inside  the  copper 

pc 

V  •  H  =  0 
Outside  the  copper 

V^H  =  0 


V  X  H  = 


77  47rJ 


V  •  H  =  0 


V  •  H  =  0 


The  current  normal  to  any  copper- air  surface  must  be  zero. 

The  eddy  currents  in  the  sense  lines  produce  a  negligible  field,  since  the  sense  lines  are  so 
small.  Therefore,  we  ean  approximate  our  structure  by  that  shown  in  Fig.  3.  The  eddy  currents 


Fig.  3.  Capper  structure  far  eddy-current  calculations. 
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in  the  bit  itself  are  assumed  negligible,  since  it  is  extremely  thin.  The  easy  axis  is  in  the  y- 
direction.  The  striplines  extend  to  plus  and  minus  infinity  in  the  y-direction.  As  a  further 
simplification,  we  look  only  for  those  modes  which  are  symmetrical  about  the  z  -  0-plane, 
which  amounts  to  assuming  that  the  film  is  in  the  exact  center  of  the  stripline  structure. 

The  problem,  as  stated,  is  fairly  difficult.  First,  we  will  solve  the  simpler  problem  in 
which  the  striplines  are  replaced  by  infinite  copper  sheets,  and  W  goes  to  infinity.  We  can 
solve  the  partial  differential  equations  involved  by  separation  of  variables.  We  assume  a  solu¬ 
tion  of  the  proper  form  for  each  of  the  three  regions  involved  (copper,  air  between  lines,  air 
outside  of  lines)  and  apply  boundary  conditions.  Region  one  is  the  volume  outside  the  lines. 
Assume  that 

exp[— a'z]  sin  (by)  cos  (gx) 


In  order  to  satisfy  the  divergence  condition  on  H  (that  it  be  zero), 
Hx  =  H^x  exp[-a'z]  sin  (by)  sin  (gx) 

Hy  =  exp[— a'z]  cos  (by)  cos  (gx) 

Applying  the  condition  that  the  divergence  equal  zero, 

-a'lF  -bH’  +  glF  =0 
oz  oy  ^  ox 


Since  the  Laplacian  of  II  must  also  equal  zero, 


In  region  two  (the  volume  inside  the  copper), 

ll^  =  sin(a’'z  +  n)  sin  (by)  cos(gx) 

Ilx  -  cos(a'’z  +  n)  sin  (by)  sin(gx) 


=  ^^oy  ^ 

a”H'’  +  glF'  -  bll”  =  0 

oz  ^  ox  oy 

,,2  ,  ,  2  ^  2  47rs 

a"  +  b  +  g  = - 2 

pc 

In  region  three  (the  volume  between  the  striplines), 

II^  =  sinh(a*’'z)  sin  (by)  cos(gx) 

Ilx  "  ^^ox  cosh  (a’ "z)  sin  (by)  sin(gx) 

=  ^^oy  cosh  (a' "z)  cos  (by)  sin(gx)  , 

,1.2  k2  ,  2  ,2 

a'”  =  b  +  g  =  a'  , 

a'lF'’  +  gIF”  -  IF”  -  0 
oz  ^  ox  oy 


The  z -component  of  the  current  must  equal  zero  at  z  =  h  and  z  =  h  +  jP  . 
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Either 


an 


an. 


X 


z'"  ax  ay  ^ 

0  -  (H”  g  +  bH’ '  )  cos  (a*  'z  +  n)  cos  (by)  sin  (gx) 


cos  (a*  'h  +  n)  =  0 


cos  [a'  ’(h  +  i  )  +  n]  =  0 


or 


copper-air  boundaries. 

At  the  boundaries,  H  must  be  continuous.  This  last  condition  gives  us  a  sot  of  simultane¬ 
ous  equations  which  we  can  solve  (see  Appendix  A).  Deriving  an  equation  for  the  relation  be¬ 
tween  a'  and  a”,  we  find 


a 'a'*  [tanh  (a’h)  +  1]  =  [a"^  tanh  (a'h)  —  a'^]  tan  (a”  f) 


(3) 


This  equation  can  be  solved  numerically.  Note  that  a”  is  a  multivalued  function  of  a'.  We  may 
now  plot  s  as  a  function  of  a’  and  a”.  A  set  of  normalized  curves  is  presented  in  Eig.4.  l-'or 
given  values  of  b  and  g,  the  mode  with  the  lowest  value  of  s  is  called  the  fundamental  mode. 
Note  that  the  first  higher  order  mode  has  a  much  greater  value  of  s  than  the  fundamental  mode. 
Since  s  is  the  inverse  of  the  time  constant,  the  higher  order  modes  (and  modes  with  large  values 
of  b  and  g)  die  out  very  quickly.  Therefore,  after  a  short  time  the  characteristic  lime  <’onslant 
of  decay  will  be  the  value  for  the  fundamental  mode  with  a'  equal  to  zero.  Therefore, 


l/h  a' '  tan  (a"  H 


if  a" 6.  is  small.  'Hiis  turns  out  to  be  true.  Therefore,  to  a  good  approximation. 


('onsider  a  sudden  change  of  bit  magnetization  in  tlie  y-direction  at  time  zero.  At  time  Z(‘ro- 
plus,  the  field  in  the  copper  has  not  changed,  so  we  must  write  the  difference  between  the  static’ 
bit  field  in  the  copper  before  and  after  the  change  (the  transient  field)  as  a  superposition  of  the 
normal  modes  which  we  have  found.  Since  we  know  the  time  dependence  of  each  normal  mode, 
we  can  find  the  time  dependence  of  the  eddy  field  as  it  decays  toward  zero.  In  particular,  we 
are  very  interested  in  the  magnitudes  of  the  lower  modes,  since  the  higher  ones  quickly  die  out. 
The  lower  modes  represent  slow  variations  of  the  fields  with  the  space  variables,  b'igure  S  is 
a  sketch  of  the  fundamental  and  first  higher  order  modes  for  small  values  of  b  and  g.  I’or  th(' 
static  bit  field,  the  variation  of  between  z  h  and  z  h  ^  f  is  not  particularly  great  if  li  is 
greater  than  f  .  For  the  nonfundamental  modes,  the  variation  of  with  z  is  extreme.  W'e 
therefore  expect  the  coefficients  of  the  nonfundamenlal  modes  to  be  small  in  our  linear  super¬ 
position.  The  variation  of  with  x  and  y,  however,  is  much  moi’e  extreme,  since  most  of 
the  static  bit  field  penetrates  the  copper  in  a  small  region  around  the  poles  of  the  l)it.  There¬ 
fore,  the  coefficients  of  modes  with  large  values  of  b  and  g  are  large.  This  tells  us  that  mu(  h 
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Fig.  4. 


Graph  of  narmalized  inverse  time  constant  (Sn)  vs  normal  ized  wavelength. 
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of  the  flux  participates  in  a  fast  decay  process,  and  only  a  small  j^ortion  of  the  flux  goes  into  tlie 
slowly  decaying  modes. 

This  makes  much  physical  sense.  The  high  concentration  of  the  static  bit  field  near  the 

poles  of  the  bit  requires  a  high  concentration  of  current  near  the  poles  to  keep  the  field  i^enetrat- 

ing  the  copper  at  time  zero-plus  the  same  as  it  was  at  time  zero-minus.  This  j^roduces  very 
2 

high  J  p  losses,  a  rapid  rate  of  dissipation  of  energy,  and  therefore  a  rapid  decay.  In  tlie  modes 
with  small  values  of  a’,  the  current  distribution  is  much  less  concentrated  and,  therefore,  the 
energy  dissipation  and  decay  rate  are  much  smaller. 

An  attempt  was  made  to  find  the  coefficients  of  the  modes  by  hand  calculation,  it  was  con¬ 
cluded  that  an  adequate  calculation  could  only  be  done  by  comi:)uter,  since  the  characteristic 
"spectrum"  for  a  reasonable  magnetic  thin  film  configuration  is  very  comjilex. 

A  similar  calculation  can  be  performed  for  a  single  sheet  of  copper.  We  find  tliat 


a 


,2 


a”  tan(a*'i/2) 


47rs 

2 

pc 


(5) 

(6) 


For  a’  smaller  than  l/h,  which  is  the  case  for  the  slowly  decaying  modes,  the  single-sheet  time 
constant  is  much  longer  than  the  double-sheet  time  constant.  This  again  makes  })hysical  sense, 
since  the  field  in  the  single-sheet  case  is  much  more  extensive  than  in  the  double-sheet  case. 
Therefore,  there  is  more  field  energy  stored  in  the  single-sheet  H  field  tlian  in  tlie  double -slu'et 
II  field  for  the  same  current  densities  and  energy  dissipation  rates. 

As  the  width  of  the  stripline  apj^roaches  the  width  of  the  bit,  the  validity  of  tlie  infinite  sheet 
tlu^ory  becomes  questionable.  However,  the  theory  is  still  not  comjjletely  unreasonable,  since 
the  aspect  ratio  between  the  line  widths  and  their  S})acing  is  still  quite  large.  It  would  be  very 
desirable,  however,  to  have  a  theory  that  exj^licitly  took  the  width  of  the  stri})lines  into  account. 

To  make  the  problem  tractable,  we  assume  that  the  striplines  are  thin  sheets  and  tliat  tlie 
current  is  uniform  across  them.  This  means  that  only  the  time  derivative  of  H.^  is  effective  in 
producing  a  current.  Because  of  this  assumption,  we  only  find  modes  similar  to  the  fundamental 
modes  of  the  infinite  sheet  theory.  Once  again,  when  we  have  found  the  normal  modes,  we  will 
write  as  a  linear  superposition.  Many  techniques  were  tried  in  order  to  solve  the  problem. 
The  one  finally  chosen  is  fairly  straightforward,  though  laborious. 

We  can  find  the  value  of  11^  by  integration  over  tlie  currents  of  the  sheet.  We  consider  J  as 
a  surface  current,  and  p  as  a  surface  resistivity  (bulk  resistivity/thickness) ; 


II 

z 


1) 

z 


JyAy  -  y')  -  Jy(x  -  x') 
[(y  -  y')^  +  (x  -  x')^  +  (z  - 


dx'  dy' 


(V  X  J)  •  1^  (at  the  planes) 


(7) 


The  integration  must  be  performed  over  both  the  upiier  and  lower  {ilanes.  We  now  assume  a 
form  for  the  currents  in  tlie  plane.  We  assume  that  the  current  in  the  lower  plane  is  the  negative 
of  the  current  in  the  Ujiper  jilane,  by  symmetry.  Since  11  must  be  a  symmetric  function  of  x, 
the  current  in  the  y -direction  must  be  antisymmetric  in  x.  The  x-comiionent  of  the  current  must 
be  zero  at  x  -  ±W/2.  We  chose  to  represent  J  by  a  trigonometric  series,  so  J  goes  as  sin(Nx). 

^  J 

Sin(n7rx)  forms  a  comjilete  set  from  zero  to  one;  sin(2n  -f-  1)  ttx  is  symmetric  about  2,  while 
sin(2nx7r)  is  antisymmetric  about  I-  We  move  the  origin  of  the  set  to  \  for  convenience. 
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An  assumed  J.  which  satisfies  the  symmetry  requirements  in  the  upper  sheet,  has  its 
divergence  equal  to  zero,  and  has  a  J  , -component  equal  to  zero  at  x  -  ±W/Z  is 


where 


and 


J  L  sin(Nbx/g)  sin  (by)  +  l^(g/N)  cos(Nbx/g)  cos  (by)] 

n  0 

M  (Zn  ■+  1)  TT 

g  bW 

OO 

I.^(vx  J)  1^*  cos(Nbx/g)  sin  (by) 

n-  0 


(8) 


(9) 


Note  that  we  have  not  summed  over  b.  This  assumes  that  a  normal  mode  has  a  sinusoidal  varia¬ 
tion  in  tiie  y -direction.  To  show  that  this  is  true,  consider  the  integral  for  the  field  in  the  upper 
half  plane  produced  by  the  currents  in  the  lower  half  plane  (note,  m  -  Zh), 

foo  p  W/Z  [sin(Nbx’/g)  sin  (by')  I  (x  —  x')  —  ^  cos(-^^— )  cos  (by') 


n-O 


w/z 


dx 


[(y  -y')^  +  (X  -  x')'^  I  n/]^^ 


(10) 


("on skier  tlie  y -dependence  of  the  first  term  of  the  integrand.  Make  the  transformation 
y  -  y’  -  V ,  and  dy’  dr.  Thus, 


f-oo 


sin  (i)y ')  dy  ’ 


[(y  -  y  )  ^  (x  -  x’)  +  m  ]  ^ 


-  sin (by) 


— OO 


cos(br)  dr 


r  Z  ,  ,  ,,Z  ,  Z,3/Z 

[r  +  (x  -  X ')  I  m  ]  ' 


(11) 


and 


r" 


(y  -  y')  <-'os  (by')  dy' 

r/  1  /  .72  ,  2TJ72' 

I  (y  "■  >  )  M  X  ~  X ' )  +  m  ]  ' 


*  4-00 


sin  (by) 


r  sin  (br)  dr 


'"-«>  +  (x  -  x')'^  f 


(IZ) 


By  a  change  of  v'ar tables,  we  may  write  tiie  above  expires sions  so  tiiat  they  involve  integrals 
with  only  one  variable  parameter: 

A  \[ m^  r  (x  —  x’ )‘^  , 


and 

B  r/A 

Kquations  (11)  and  (IZ)  become,  respectively. 


sin  (by)  cos  (bAIA  dL  sin  (by)  ZCo(bA) 

^  ,,  ,  ,  ZZ/Z  “  ,Z 


(13) 


(1  ^  IZ) 


and 


sin  (by 


^  I,  sin(bXL)  dL 

(1  + 


sin  (by)  ZSo(bA) 

A - 


(H) 
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where 


L  sin  (bXL)  dL 


cos  (bXL)  dL 


Note  that  the  derivative  of  Co(x)  with  respect  to  x  is  — So(x). 

Set  (7)  equal  to  (V  X  J)  *  l^pc/s.  We  now  multiply  both  sides  of  the  expression  by  eos(Qbx/g) 
and  integrate  from  —  W/2  to  +W/2.  Since  the  cosines  are  orthogonal  in  the  interval  from  —  W/Z 
to  +W/2,  we  find  that,  after  dividing  out  the  sin  (by)  dependence  from  both  sides  and  some  manip¬ 
ulating  (for  full  details  see  Appendix  B),  the  following  expression  must  hold  for  all  values  of  q: 


(1^) 


n-0 


where 


and 


Q  =  (2q  +  1)  TT 


lower  sheets  to  the  field  at  the  upper  sheet,  and  (b)  arises  from  the  eontribution  of  the  eur- 

r  r 

rents  in  the  upper  sheet  to  the  field  at  the  upper  sheet. 


The  expression  is  described  by  one  expression  when  n  equals  q,  and  by  another 

expression  when  n  is  not  equal  to  q: 


+  2Q  sin  k)  sin  ^  k)  De) 


(16) 


(17) 


where 


Kk)  -  -  cos  (■^)  +  sin  (-^) 


k  =  (X  -  X')  b 


Ds(k,  b)  -  b\^So(k)  -  kXSo(bX) 


nc(k,  b)  -  ^  Co(bX)  -bX^Co(k) 
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For  any  value  of  b,  the  set  of  simultaneous  equations  deseribed  by  (15),  (16),  and  (17) 
have  nontrivial  solutions  only  for  eertain  values  of  s.  Corresponding  to  eaeh  one  of  these 
values  of  s  and  b,  there  is  a  normal  mode  of  deeay.  We  ean  solve  for  the  values  of  s  which 

4 

give  nontrivial  solutions  if  we  put  (15)  in  matrix  form,  i.e., 

M  A  =  0 
qn  n 


where 


M  =H  (b)-5  [1  +  (g  )^] 

qn  nq'  '  qn  Zs  ^  Q  ^ 


(18) 


qn 


^  n 


1  q  -  n 


There  is  a  nontrivial  solution  for  the  set  of  equations  only  where  the  determinant  of  is 
equal  to  zero.  The  applieation  of  this  eondition  will  give  us  a  polynomial  in  s  to  solve.  For 
eaeh  root,  we  ean  solve  for  the  eoeffieients  whieh  deseribe  the  normal  mode. 

In  order  to  make  use  of  the  results  obtained,  we  must  evaluate  Co(x).  By  differentiation, 
So(x)  may  be  found.  Therefore,  a  numerieal  integration  was  performed  and  a  few  points  plotted. 
It  was  found  that  an  exeellent  fit  eould  be  obtained  with  a  funetion  eonsisting  of  the  sum  of  two 
exponentials  (see  Fig.  6).  For  a  more  thorough  evaluation,  it  is  possible  to  write  Co(x)  as  a  sum 

5 

of  I.egendre  funetions.  The  Legendre  funetions  form  a  eomplete,  orthonormal  set,  from  zero 
to  infinity,  with  the  general  expression 

i=n 

S  ^ni  ^ 

i-0 
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The  coefficients  are  tabulated  by  Lee  and  others.  If  we  wish  to  express  Co{x)  as  a 
linear  superposition  of  with  coefficients  we  find  by  integration,  i.e,. 


Co{x)  V^{x)  dx 


ry  CO  r%  CO 

=  J  y  - 


D  =  r 

Jo 

Doing  the  x  integral  first, 

p  oo  w. 

\  exp[— w.x]  cos  (xy)  dx  =  — ^ - 2 

^  r\  ^  yir  4-  \r 


r  C^.  cxp[— px(i  +  '^)]  COS  (xy)  dy  dx 


w.  +  y 
1 


(1  + 


(ZO) 


where 


Equation  (20)  now  becomes 


W.  =  p(i  +  y) 


D  ■=  ^  C 
n  ^  ni 

i-O 


o  +  y^)  (1  +  y“) 


w.  dy 

— - TTTT  =  7  C  .F(w.) 

,  ,  2,3/2  ^  ni  '  1' 

1  +  y 


i=0 


(21) 


The  integral  can  be  evaluated  in  closed  form  upon  the  substitution: 

u 

y  "  -  - 

n/i  -u^ 


w.  dy 


o  (Wj^  +  y^)  (1  + 


y  1  w.  (1  —  u  )  du 


o  w.  +  (  1  —  w.  )  u 
1  1 


F(w.) 


w. 

1 


2  , 
w.  —  1 
1 


2  , 
w.  —  1 
1 


^  ^  In  -  1)-  In  (w.^  -  Jw.^  -  1) 

2w.  -  1 


w  >  1 


arctan 


1  - 


w. 

1 


f,  2 

V.  / 1  —  w. 

1  V  1 


w  *  1 


(22) 


Knowing  F(w^)  in  closed  form,  and  knowing  we  can  find  the  coefficients  D^;  p  is  chosen 

to  give  rapid  convergence.  This  technique  would  be  quite  suitable  for  machine  computation.  It 
is  inconvenient  for  hand  computation  because  the  functions  involve  the  small  differences  of 
large  quantities. 

A  hand  calculation  was  done,  basically,  to  see  if  the  results  obtained  were  reasonable. 

Co(x)  was  approximated  by 

Co(x)  =1.35  exp  f-0.804x]  -  0.35  exp  [-3.10x] 

This  expression  has  the  proper  value  and  slope  at  x  =  0  and  the  correct  asymptotic  behavior  for 
large  x.  Wc  choose  W  to  be  50  mils  and  h  to  be  7  mils.  The  wavelength  of  the  sinusoidal 
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variation  in  the  y-direction  is  chosen  to  be  100  mils.  The  bits  wc  use  arc  roughly  50  mils  long; 
choosing  this  value  of  b  puts  the  adjacent  positive  and  negative  peaks  of  the  sinusoid  over  the 
north  and  south  poles  of  the  bit.  We  only  consider  the  bottom  two  modes.  Using  an  interval  of 
0.1  g,  we  do  a  numerical  integration  to  calculate  the  values  of  the  matrix  components.  The  re¬ 
sultant  matrix  equation  is 


0 

0 


0.503  -  t{0.151) 
+0.103 


+0.034 

0.557  -  t(0.251) 


t 


2 


47rf  s 


1^0 


time  constant 


47rf  t 
Z 


Pb^ 


The  bulk  resistivity  is  p^,  the  surface  resistivity  is  p,  and  the  thickness  of  the  copper  is  I. 

The  two  solutions  for  t  are  3.41  and  Z.14.  I"or  convenience,  t  is  measured  in  mils.  Cor¬ 
responding  to  the  3.41  solution,  we  find  over  to  equal  +0.34  5.  Corresponding  to  the  2.14 
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Fig.  7.  Normal  modes  of  thin,  finite 
stripline  theory. 


solution,  A^  over  A^  is  equal  to  —5.3.  Assuming  the  thickness  and  resistivity  to  be  the  same 
as  in  our  structure,  we  find  the  time  constants  to  be  172  and  104  nsec.  Plots  of  the  x  variation 
of  at  the  upper  plane  are  shown  in  Fig.  7.  They  are  quite  different  from  the  modes  of  decay 

one  finds  for  the  infinite  sheet  case. 

For  b  =  0,  the  equations  can  be  put  in  a  much  simpler  form: 
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(23) 


OO 


A  =  y  A  H  (0) 
q  ^  n  nq 


n  =  0 


H  (0)  =  8m^W 
nq' 


(24) 


o 


z  z  z 

X  A  (N  -  Q^) 


g  dx{W  -  x)  sin  (■^) 


(25) 


x\' 


2 


The  ratio  between  the  widths  of  the  striplines  and  their  spacing  is  large  enough  (for  all  our 
configurations)  so  that  the  off-diagonal  matrix  elements  are  small  enough  to  make  a  negligible 
change  in  the  time  constant.  For  example,  if  we  merely  set  the  diagonal  elements  of  our  matrix 
equal  to  zero,  we  would  have  obtained  values  of  3.34  and  2.22  for  t.  This  compares  very  closely 
with  the  more  accurate  values  of  3.41  and  2.14. 

Taking  only  the  diagonal  matrix  elements, 


T  "  time  constant  - 

q 


(26) 


It  is  important  to  note  that  the  modes  we  find  for  the  same  value  of  b  are  not  orthogonal  to 
each  other. 

IV.  EXPERIMENTAL  APPARATUS 

In  the  experiments  which  were  performed,  driving  currents  were  generally  on  the  order  of 
amperes,  while  the  sense  line  output  was  on  the  order  of  a  fraction  of  a  millivolt.  The  rise  and 
fall  times  of  the  drive  pulse  were  quite  short.  As  a  result,  the  noise  problem  was  exceedingly 
severe.  With  careful  alignment  and  shielding,  however,  it  was  possible  to  increase  the  signal- 
to-noise  ratio  to  three -to-one. 

The  configuration  used  is  shown  in  Fig.  8.  The  bit  line  and  the  top  word  line  are  cut  away 
to  show  the  sense  line  and  the  bit.  A  side  view  of  the  sandwich  is  shown  in  Fig.  9.  The  dummy 
bit  line  was  not  connected  electrically  to  any  power  source.  It  was  used  to  obtain  an  estimate 
of  the  effect  of  the  eddy  currents  flowing  in  the  bit  line  upon  the  bit.  It  changed  the  measured 
eddy  field  by  only  2  percent. 

Three  different  word  lines  were  used.  The  7  5-  and  150 -mil-wide  lines  were  both  made 
from  the  same  material.  The  thickness,  as  measured  by  a  micrometer,  was  2  mils.  The  more 
important  parameter  for  our  work,  however,  was  the  surface  resistivity.  This  was  found  by 
measuring  the  resistance  of  a  length  of  the  stripline,  and  multiplying  it  by  the  ratio  of  width  to 
length.  The  BO-mil-wide  stripline  was  made  of  different  material.  Its  thickness  also  measured 
2  mils,  but  its  surface  conductivity  was  somewhat  lower. 

The  sense  line  is  different  from  the  one  shown  in  Fig.  1.  This  type  of  sense  line  is  known  as 
a  planar  sense  line.  It  reduces  the  number  of  layers  required  between  the  word  lines,  so  that 
they  may  be  brought  closer  together.  This  type  of  sense  line  also  decreases  the  coupling  be¬ 
tween  the  bit  line  and  the  sense  line.  With  the  structure  shown  in  Fig.  1,  a  very  large  signal 
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Fig. 8.  Experimentol  stripl ine  configurotion , 
top  view. 


Fig,  9,  Experimentol  stripline  configurotion, 
side  view. 


Fig.  10.  Associoted  electronics. 
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is  coupled  into  the  sense  line  during  a  write  operation,  and  the  sense  amp^lifier  may  not  be  used 
again  until  it  has  recovered  from  the  effects  of  the  large  pulse.  The  total  thickness  of  the  sense 
line  and  backing  was  found  to  be  4  mils. 

The  thin  film  plane  used  contained  square  bits  50  mils  on  a  side.  The  inter-bit  spacing  was 
also  50  mils.  The  plane  was  originally  7  mils  thick,  giving  a  total  word-line -to-word -line  spac¬ 
ing  of  11.25  mils.  For  the  last  part  of  the  experiment,  the  glass  plane  was  etched  until  it  was 

only  3  mils  thick,  giving  a  total  spacing  of  7.2  5  mils. 

The  driving  and  sensing  circuitry  is  shown  in  Fig.  10.  The  sense  lines  are  connected  to 

the  sense  amplifier  with  double -wire,  shielded  cable.  The  sense  amplifier  is  a  balanced  differ¬ 

ence  amplifier  with  a  gain  of  100.  The  output  of  the  sense  amplifier  is  displayed  on  a  Type  535 
Tektronix  oscilloscope  which  has  a  53B  plug-in  unit.  The  rise  time  of  the  sense -amplifier- 
oscilloscope  system  was  about  20  nsec.  The  sync  amplifier,  negative-bit  driver,  and  slow  word 
pulse  generator  used  were  all  Lincoln  Laboratory  Model  V  standard  pulse  generators.^  They 
are  capable  of  driving  a  maximum  of  about  500  ma  into  a  50 -ohm  load,  with  rise  and  fall  times 
on  the  order  of  150  nsec.  The  positive -bit  pulse  driver  was  a  Lincoln  Laboratory  Model  VI 
pulse  generator.  The  two  types  are  identical,  except  that  the  Model  V  delivers  a  negative  pulse, 
and  Model  VI  delivers  a  positive  pulse. 

The  fast  pulse  generator  is  a  Kese  Type  1051  millimicrosecond  pulse  generator.  At  maxi¬ 
mum  output,  it  can  drive  a  50-ohm  load  with  a  4-amp  pulse  having  rise  and  fall  times  of  15  and 
2  5nscc,  respectively.  For  a  2 -amp  setting,  the  rise  and  fall  times  are  10  and  18  nsec,  respec¬ 
tively.  The  pulse  widths  are  selected  by  a  switch.  The  available  pulse  widths  arc  10,  20,  50, 
and  lOOnscc. 

The  timer  allows  us  to  operate  the  equipment  in  a  very  flexible  manner.  A  single  timer 
cycle  contains  eight  intervals.  Any  combination  of  pulse  drivers  and  scopes  may  be  turned  on 
during  each  interval.  The  basic  cycle  frequency  is  limited  to  about  200  cps,  because  of  the  long 
recovery  time  of  the  fast  pulse  generator.  An  interval  stretch  option  is  available  which  allows 
us  to  repeat  one  interval  several  times  before  continuing  with  the  cycle. 

Some  of  the  data  on  the  film  were  taken  with  a  Laboratory-built  H-II  looper.  This  device 
measures  the  total  magnetization  of  the  film  as  a  function  of  the  driving  field,  and  |)resents  the 
result  on  an  oscilloscope.  The  axes  of  the  driving  field  and  measured  magnetization  may  be 
either  parallel  or  perpendicular  to  each  other. 

Rather  than  measure  k,  we  normally  measure  a  quantity  H^,  which  is  the  field  required  to 

completely  rotate  the  magnetization  into  the  hard  direction.  For  a  field  less  than  11^  (if  we 

apply  a  field  in  the  hard  direction  only),  equating  M  x  II  to  the  anisotropy  torque  gives  ull^  cos  A 

2k  sin  A  cos  A.  This  means  that  where  A  first  reaches  90®,  =  2k/u.  For  the  plane  used, 

H,  4.3  oe. 
k 

To  take  the  measurement,  the  field  is  applied  in  the  hard  direction,  and  the  magnetization 
is  measured  in  the  hard  direction.  The  saturation  of  the  flux  in  the  hard  direction  is  spread  out 
due  to  dispersion  so  that,  usually,  is  taken  to  be  the  value  of  the  field  at  which  a  line  drawn 
tangent  to  the  B  vs  H  curve  at  the  origin  intersects  the  saturation  value  of  13. 

Dispersion  may  also  be  measured  with  the  I3-H  looper.  The  plate  in  the  device  on  which  the 
film  is  mounted  can  be  rotated.  If  the  driving  field  is  roughly  parallel  to  the  easy  axis,  and  the 
sensing  coil  perpendicular  to  the  drive  coil  (see  Fig.  11),  the  magnetization  measured  by  the 
looper  is  the  difference  between  the  portion  of  the  magnetization  to  the  left  and  the  portion  of 
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Fig.  1 1.  B-H  looper  plote  mount. 


Fig.  12.  Bias  field  vs  sense  omplifier  output 
hysteresis. 
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the  magnetization  to  the  right  of  the  field  axis.  This  is  because,  under  the  influence  of  the 
field,  the  two  portions  of  the  magnetization  rotate  in  opposite  directions,  and  their  components 
along  the  sense  eoil  axis  cancel  out.  The  difference  between  the  two  angles  where  the  total 
magnetization  rotates  as  a  unit  is  the  dispersion  angle.  Normally,  the  angle  between  the  90- 
pereent  points  is  measured. 

The  magnetization  of  the  film  is  a  constant  of  the  material.  It  is  measured  by  a  ferromag¬ 
netic  resonance  technique  and  has  been  shown  to  be  10,000  oe  for  the  material  we  use.  The 
composition  of  the  material  is  83 -percent  iron  and  17 -percent  nickel. 

The  film  was  evaporated  on  optical  cover  glass  in  a  DC  field  of  23  oe.  It  was  evaporated 

o  o 

for  30  seconds  at  an  evaporation  rate  of  1000  A  per  minute,  giving  a  thickness  of  500  A.  The 
field  necessary  to  produce  extensive  wall  motion  in  the  easy  direction  is  1  oe. 

V.  EXPERIMENTAL  PROCEDURE  AND  RESULTS 

1 

At  present,  there  is  little  data  eoneerning  eddy  eurrent-NDRO  in  the  literature.  There¬ 
fore,  rather  than  take  extensive  sets  of  readings  for  a  single  system  (varying  only  one  or  two 
parameters)  it  seemed  wiser  to  concentrate  on  varying  several  of  the  parameters,  taking  less 
extensive  data  for  each  one.  Two  types  of  experiments  were  made. 

The  first  experiment  was  performed  to  demonstrate  the  possibility  of  eddy-NDRO.  We  con¬ 
tinually  pulsed  the  word  line,  monitoring  the  output  of  the  sense-amplifier  and  varying  the  DC 
bias  produced  by  the  bit  line.  If  we  started  with  a  high  positive  value  of  bias  and  gradually  de¬ 
creased  it,  the  bit  would  initially  be  set  in  the  positive  direction.  Making  the  bias  negative 
enough  will  reverse  the  direction  of  the  bit.  The  same  procedure  ean  be  followed  going  from  a 
negative  to  a  positive  bias.  When  the  results  were  plotted,  a  hysteresis  loop  was  observed. 

This  means  that  there  is  a  range  of  operating  conditions  over  which  the  system  output  has  two 
stable  values  (depending  on  history),  and  NDRO  is  possible  (see  Fig.  12). 

The  signals  for  large  bias  levels,  for  which  NDRO  is  impossible,  correspond  to  complete, 
coherent  rotation  of  flux.  For  the  larger  loops,  there  is  very  little  decrease  in  output  in  the 
eddy-NDRO  mode.  The  smaller  loop  displays  a  greater  signal  loss.  (We  are  probably  approach¬ 
ing  the  limit  at  which  NDRO  ean  be  obtained  in  this  ease.) 

When  the  bias  was  reduced  to  zero,  the  signals  appeared  perfectly  stable,  with  no  measur¬ 
able  change  after  several  million  read  operations.  NDRO  at  zero  bias  eould  be  obtained  with  all 
the  configurations  that  were  tried.  Although  no  detailed  data  were  taken,  the  effects  of  bit  re¬ 
generation  were  noted.  If,  for  example,  a  bit  set  in  the  positive  direction  was  given  a  negative 
bias  of  sufficient  value  to  cause  the  signal  output  to  decrease  appreciably  (due  to  break-up),  the 
output  signal  would  increase  again  as  the  eddy  current  field  regenerated  the  bit  —  upon  reduction 
of  the  magnitude  of  the  negative  bias. 

The  second  experiment  was  done  to  obtain  more  detailed  information  about  the  switching 
process.  We  attempted  to  measure  the  eddy  field  as  a  function  of  time.  This  is  a  rather  diffi¬ 
cult  measurement  since  the  fields  are  small,  decay  rapidly,  and  coexist  with  the  word  drive 
field  (which  is  much  larger  in  magnitude).  The  method  finally  chosen  is  rather  indirect.  Con¬ 
sider  a  bit  which  is  set  in  the  positive  direction.  With  a  small  DC  field  in  the  negative  direction, 
a  word  pulse  is  applied.  At  the  end  of  the  word  pulse,  the  effective  field  in  the  easy  direction  is 
the  difference  between  the  eddy  and  the  bias  fields.  If  the  eddy  field  is  larger  than  the  bias  field, 
the  magnetization  will  fall  back  in  the  positive  direction  and  the  sense  signal  on  the  trailing  edge 
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Results 


System  Porometers 


Time  constont  76  nsec 

90-percent  transition 
width  0.  38  oe 

Normolized  coefficient  0.  143 


Spacing 
Line  width 
Driving  field 


7.  25  mils 
50  mils 
17.8  oe 


Fig.  13.  Bias  field  vs  sense  amplifier  output,  eddy  field  meosurement. 


Fig.  14.  Normalized  time  constant  for  thin, 
finite  striplines 


Fig.  15.  Inverse  normalized  time  canstant  for  thick, 
infinite  sheet. 
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of  the  word  pulse  will  have  a  polarity  opposite  to  that  of  the  signal  on  the  leading  edge  of  the 
word  pulse.  If  the  bias  field  is  larger,  the  polarities  will  be  the  same.  If  the  bit  were  set  in 
the  negative  direction,  a  similar  transition  would  occur  for  positive  values  of  bias.  Measuring 
the  difference  between  the  two  transitions  gives  twice  the  effective  eddy  field  (see  Fig.  13). 

It  would  have  been  better  to  apply  the  field  in  the  easy  direction  only  for  the  short  time  that 
the  word  drive  was  on,  rather  than  apply  it  in  a  DC  manner.  This  is  because  some  wall  motion 
may  occur  during  the  time  between  the  writing  and  reading  of  the  pulse,  due  to  the  influence  of 
the  DC  bias  field.  The  coupling  between  the  bit  and  sense  lines  is  large  enough,  however,  so 
that  the  sense  amplifier  does  not  recover  for  a  considerable  period  of  time  after  the  application 
of  the  easy  direction  field,  if  it  is  done  in  a  pulsed  manner.  The  effects  of  wall  motion  would 
be  indicated  by  a  decrease  of  the  signal  obtained  on  the  leading  edge  of  the  word  pulse  with  an 
increase  in  bias.  This  effect  was  not  noted  for  fields  on  the  order  of  those  needed  to  produce 
the  transitions. 

The  measurement  suffers  from  several  inaccuracies.  The  poor  signal-to-noise  ratio  made 
the  reading  of  the  output  signal  value  difficult.  Also,  the  dispersion  of  the  bit  spread  out  the 
transition.  The  noise  could  be  lowered  by  reducing  the  word  drive.  Reducing  the  word  drive, 
however,  makes  the  assumption  that  the  bit  magnetization  is  completely  switched  into  the  hard 
direction  on  the  leading  edge  of  the  word  drive  pulse  less  valid.  The  value  chosen  was  a  com¬ 
promise.  From  Fig.  12,  it  can  be  seen  that  there  is  little  change  in  the  hysteresis  loops  when 
the  driving  field  is  increased  from  17.6  to  29.8  oe.  The  shape  of  the  loops  is  greatly  affected  by 
the  size  of  the  eddy  field.  This  verifies  the  fact  that  the  easy-to-hard-direction  switching  time 
is  not  important  for  fields  much  larger  than  Hj^,  which  is  4.3  oe.  For  drive  fields  of  the  order 
of  FIj^,  the  output  signal  does  decrease. 

There  is  some  error  due  to  the  neglect  of  the  interaction  between  neighboring  bits.  This 
error  amounts  to  only  0.03  oe,  or  less,  and  would  be  a  constant  for  any  experiment,  i.e.,  it 
merely  shifts  the  zero. 

The  rise  time  of  the  driving  pulse  is  not  zero.  However,  our  times  are  the  differences  be¬ 
tween  the  50-  and  100-nsec  readings.  The  correction  for  the  finite  rise  and  fall  times  will  in¬ 
crease  both  values,  keeping  the  difference  roughly  constant. 

-4 

The  measured  values  of  p  were  found  to  be  6.5  X  10  ohm  per  square  for  the  50-mil  lines, 
-4 

and  5.37  x  10  ohm  per  square  for  the  150-  and  75-mil  lines.  Using  these  values,  we  calculate 
the  theoretical  results  shown  in  Table  1.  The  theoretical  values  for  the  thin,  finite -width  strip¬ 
line  approximation  b  -  0  were  found  from  (25)  and  (26)  using  Fig.  14  (which  is  a  normalized  plot 
of  these  equations).  The  theoretical  values  for  the  thick,  infinite -width  sheet  approximation 
a'  =  0  are  found  from  (3)  using  Fig.  15. 

The  experimental  results  are  shown  in  Table  II.  The  estimated  accuracy  of  the  individual 
eddy  field  measurements  is  ±10  percent,  and  of  the  time  constant,  ±30  percent.  It  is  expected 
that  the  actual  time  constant  of  the  slowest  decaying  mode  will  be  greater  than  the  value  calcu¬ 
lated  by  the  thin,  finite -width  sheet  approximation  and  smaller  than  the  value  calculated  by  the 
infinite -width  sheet  approximation.  The  measured  time  constant  would  be  less  than  the  actual 
time  constant  of  the  slowest  decaying  mode  because  of  the  presence  of  higher  order  modes. 

This  is  seen  to  be  the  case.  The  thin,  finite -width  sheet  approximation  seems  to  be  the  better 
approximation  to  the  actual  behavior  of  the  system. 
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TABLE  1 

THEORETICAL  VALUES 

Thin,  Infinite-Width  Sheet  [Eq.(4)] 

Line  Width 
(mils) 

Separation 

b 

(mil-1) 

7.25  mils 

11.25  mils 

14  mils 

50 

178 

272 

75 

217 

336 

150 

217 

336 

Thin,  Finite-Width  Striplines  [  Eqs 

.(25),  (26)1 

50 

120 

160 

182 

0 

75 

165 

225 

0 

150 

191 

272 

0 

50 

167 

0.0628 

Thick,  Infinite-Width  Sheet,  o'  =  0  [Eq.(3)] 

50 

212 

305 

75 

259 

376 

150 

259 

376 

TABLE  11 

EXPERIMENTAL  RESULTS 

Drive 

Width 

Spacing 

Time  Constont 

Normolized 

(oe) 

(mils) 

(mils) 

(nsec) 

Coefficient 

13.5 

150 

7.25 

198 

0.202 

13.5 

75 

7.25 

120 

0.208 

17.8 

50 

7.25 

76 

0.143 

17.8 

50 

11.25 

133 

0.099 
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Along  with  the  time  constant,  the  coefficient  of  the  exponential  decay  is  required.  It  would 
be  possible  to  find  this  using  the  theory  that  has  been  presented  by  doing  the  detailed  linear 
superposition  of  normal  modes.  This  involves  a  very  extensive  calculation  which  has  not  yet 
been  made.  Instead,  let  us  obtain  an  upper  bound  on  the  coefficient.  Assume  that,  upon  switch¬ 
ing  the  bit  magnetization  into  the  hard  direction,  all  of  the  bit  flux  is  trapped  in  the  copper  by 
the  eddy  currents,  so  that  the  initial  eddy  flux  is  the  total  flux  of  the  bit.  Also  assume  that  the 
entire  flux  pattern  decays  with  the  time  constant  of  the  lowest  order  mode  (and  that  no  fast  decay 
process  takes  place).  Finally,  assume  that  there  is  no  fringing,  i.e.,  that  the  field  is  uniform 
across  the  separation  of  the  word  lines  and  does  not  extend  in  the  x-  and  y-directions  beyond  the 
borders  of  the  bit.  In  this  case  the  coefficient  of  the  exponential  decay  will  be  ud/2h,  where  d 
is  the  thickness  of  the  bit  and  u  is  its  magnetization.  The  actual  measured  value  of  the  coeffi¬ 
cient  divided  by  ud/2h  is  given  in  the  fifth  column  of  Table  II.  As  can  be  seen,  this  is  much  less 
than  1,  due  to  the  effects  of  the  fringing  of  the  eddy  field  and  the  fact  that  much  of  the  bit  flux 
engages  in  a  rapid  decay  process.  The  variation  is,  as  expected,  decreasing  for  structures  with 
large  values  of  2h/W,  and  increasing  to  a  limit  for  those  structures  for  which  there  is  a  large 
overlap  of  the  striplines  over  the  bit. 

If  the  normalized  coefficient  were  1,  the  initial  value  of  the  normalized  eddy  field  would  be: 

Spacing  (mils)  ud/2h  (oe) 

7.25  2.72 

11.25  1.75 

For  the  structure  of  two  50 -mil  word  lines  separated  by  11.25  mils,  NDRO  is  just  possible 
at  100  nsec  —  as  is  evidenced  by  the  smallest  loop  of  Fig.  12.  The  eddy  field  at  the  end  of  the 
word  pulse  in  this  case  is  0.0818  oe.  This  is  less  than  the  dispersion  from  the  easy  axis,  which 
is  0.172  oe  at  the  90-percent  point.  This  is  not  unreasonable  since  most  of  the  magnetization  is 
quite  close  to  the  easy  axis,  and  the  90 -percent  point  is  well  out  on  the  tail  of  the  dispersion 
curve.  Also,  some  breakup  does  occur  in  this  case.  We  are  able  to  obtain  another  estimate  of 
the  dispersion  from  the  width  of  the  transition.  From  the  100 -nsec  curves  associated  with  the 
structure  with  7.25-mil  spacing  and  50 -mil-wide  lines,  we  find  that  90 -percent  dispersion  from 
the  easy  axis  occurs  for  a  field  of  0.186oe.  This  agrees  well  with  the  looper  value. 

VI.  SUGGESTIONS  FOR  FURTHER  STUDY 

The  work  reported  here  must  be  considered  as  representing  an  initial  investigation  into  a 
complex  phenomenon.  Much  more  remains  to  be  done. 

To  obtain  the  maximum  amount  of  information  from  the  theory,  the  detailed  linear  super¬ 
position  of  normal  modes  should  be  performed.  It  would  be  very  interesting  to  see  if  the  low 
value  of  the  normalized  coefficient  can  be  explained  by  the  rapid  decay  of  higher  order  modes, 
as  has  been  conjectured.  The  inclusion  of  the  higher  order  modes  should  also  bring  the  meas¬ 
ured  and  calculated  values  of  the  time  constant  into  closer  agreement.  It  would  also  be  very 
useful  if  a  simple  means  of  taking  into  account  both  finite  stripline  thickness  and  width  could  be 
found. 

When  a  better  understanding  of  the  details  of  the  decay  has  been  obtained,  we  can  substitute 
back  into  (1)  to  find  the  dynamic  switching  behavior  of  the  bit.  Some  preliminary  calculations 

3 

were  performed  assuming  a  simple  exponential  decay  for  f(p).  In  agreement  with  Smay,  it  is 
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found  that  the  eddy  field  will  produce  resistive  damping  if  the  time  constant  of  the  eddy  field  is 
much  less  than  the  rise  time  of  the  driving  field.  If  the  time  constant  of  the  eddy  field  is  much 
larger  than  the  rise  time  of  the  driving  field,  the  eddy  field  has  an  instantaneous  value  which  is 
just  proportional  to  the  negative  of  the  change  in  bit  magnetization  and  has  a  clamping  effect  on 
the  switching  of  the  bit.  During  the  short  time  it  takes  the  bit  to  switch,  the  eddy  field  still 
decays  rapidly,  due  to  the  higher  order  modes.  The  effective  decay  time  might  be  short  enough 
to  bring  us  into  a  range  where  this  decay  time  is  of  the  same  order  of  magnitude  as  the  rise  time 
of  the  driving  field,  and  neither  of  the  preceding  approximations  to  (1)  are  adequate. 

The  experimental  setup  could  also  be  improved.  With  a  more  flexible  pulse  generator  and 
a  sense  amplifier  with  a  wider  bandwidth,  it  might  be  possible  to  do  a  fairly  complete  job  of 
plotting  the  eddy  field  vs  time.  The  use  of  a  bit  with  smaller  dispersion  would  also  increase 
the  accuracy  of  the  measurement.  Finally,  a  most  interesting  experiment  would  be  to  attempt 
the  construction  of  a  prototype  eddy-NDRO  memory. 
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APPENDIX  A 


Since  is  zero  in  regions  I  and  III, 


Ho’y  _  Ho^»*y  _  _  ;b  _  Ho*»y 
g  Ho”x 


Ho’x  Ho”'x 

Substituting  (A-1)  into  the  relations  derived  from  V  •  H  =  0  gives 

2 


a”Ho”z 


a’Ho”  ’z 


a’Ho’z 


Ho”x 


Ho”  ’x 


Ho’x 


The  boundary  conditions  at  z  =  h  give 

Ho”z  sin(a”h  +  n)  =  sinh(a'h)  Ho’”z 
Ho”x  cos(a”h  +  n)  =  cosh(a'h)  Ho”'x 
Dividing  (A-5)  by  (A-6)  and  using  (A-3)  and  (A-4)  gives 

tan(a”h  ^  tan(a’h) 

At  z  =  h  +  ^  we  find 

Ho’z  exp  [—  a'(h  +  f)  ]  =  Ho”z  sin  [a”(h  +  /)  +  n ] 
Ho’x  exp  (— a’(h  +  f)]  =  Ho”x  cos  [a”(/  +  h)  +  n] 
Dividing  (A-8)  by  (A-9)  and  applying  (A-2)  and  (A-3)  gives 

a” 


tan  [a”(h  +  f)  +  h]  =  — 

Using  (A-10)  and  (A-7),  we  can  eliminate  n  and  find  a  relation  between  a”,  h,  I,  and  a’: 


(A-1) 

(A-2) 

(A-3) 

(A-4) 

(A-5) 

(A-6) 

(A-7) 

(A-8) 

(A-9) 

(A-10) 


(a'/a”)  tan  (a”/) 


tanh  (a’h)  +  1 


(a”/a’)^  tanh  (a’h)  —  1 
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APPENDIX  B 


Consider  Eq.  (10).  Apply  Eqs.  (13)  and  (14)  and  divide  out  the  sinusoidal  variation  in  y.  We 
find  (10)  becomes 


2  Vl] 


pg/2  (x-x*)  Co(b\)  sin  (  ^  cos  )  So(b\) 


dx 


n=0 


-g/2 


We  multiply  through  by  cos  (Qbx/g)  and  integrate  from  —  W/2  to  W/2.  It  is  convenient  to  do 
the  integral  from  — g/2  to  g/2,  integrating  with  respect  to  u  =  bx.  Define  u'  =  bx'.  The  expression 
becomes 


[(u  -  u')  Co(bX)  Sin(^^)  -  ^  cos  So(b\) 


cos  ( 

S 


bX 


We  now  define  u  —  u'  =  k;  du'  =  — dk.  Since  X  is  only  a  function  of  K,  we  may  do  the  integral 
with  respect  to  u  in  closed  form.  Since  changing  the  sign  of  both  u  and  u'  does  not  change  the 
value  of  the  integrand,  we  need  only  consider  positive  values  of  K  and  multiply  the  answer  by  2. 
The  expression  we  must  integrate  becomes 


;  2  pg  pg/2 

Z  Hb  j, 

n=0 


dk  \ 

o  •^-g/2+k 

2  {K  sin  (u  -  K)  ]  Co(bX)  -  cos  (u  -  K)  ]  So(bX)}  cos  (^) 


bx 

It  is  elementary  to  show  that  if  q  is  not  equal  to  n: 


r  du  sin  [ ^  (u  -  K)  ]  cos  ( ^)  =  -  -I 
J-g/2  e  8 


,  Q  +  N  N 

cos(^5L_ 

ntq 


/  N  -  Q  N 
cos  ( - —  u  —  —  K) 

- g - — - g - 

N  -  Q 


g/2 

-g/2+K 


,Q  +  N  g.  ,  7r(2n  +  2q  +  2)  ,  /  i  \ 

cos  (  - -  •  =  cos  [  — - -  ]  =  —cos  7r(n  +  q) 


cos  ( 


Q  -  N,  /N  -Q,  /  ,  . 

— )  =  cos  ( — =  cos7r(n  — q) 


If  n  “  q  is  odd,  n  +  q  is  also.  If  n  —  q  is  even,  so  is  n  +  q.  Therefore, 

,N  — Q»  ,  4vn-q  .N  +  Q* 

cos( — 5-^)  =  (-1)  ^=-cos( — =-^) 


(B-1) 


.  ,N-Q,  .  ,N  +  Q,  ^ 

sin( — =  sm( — 2“^)  =  0 

Applying  these  results,  and  computing  the  value  of  the  integral,  we  find 

^  •  ,N  ,  Qu  2g(-l)""‘’Q  .  ,N  +  Q„.  .  ,Q-N„. 

•^-g/2+K  S  '  g  _N^  2g  2g 
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In  a  similar  manner,  it  can  be  shown 


’g/2 
-g/2+K 
If  n  =  q,  we  find 


r  du  cos  (u  -  K)  ]  cos  (^)  =  ^2  '  [N  sin(^|^)  -  Q  sin  (^) 

'^-cr/2+K-  S  s  w  w 


r  du  cos  {^)  sin  (u  -  K)  ]  =  ^  S  sin  ^ 

•J-a/?A-K  g  g  g 


K 


-g/2+K 

r  du  cos  (^)  cos  (u  -  K)  ]  =  ^2'^  cos  i—^)  +  ^  sin 

'^-g/Z+K  g  g  ^  g  zw  g 

We  may  now  apply  these  results  to  (B-1).  First  we  rewrite  (B-1)  as 


y.  A  H^b) 
cb  ^  n  nq 

n=0 


(B-2) 


4)  b\  (Q  -  N  )  ‘^S 


Q  -  N 


K] 


+  (Q  sin(^^)  -  N  sin(^^)]  So(b\)} 


g 


H 


K(K  -  g)  -  /v,,\  •  /QK, 

J^(b)  =  4  £^K  <  J)  -  ^ 


gb\ 


So(b\)  I(K) 


I(K)  =  cos(^)  +  ^  sin(^)  . 

There  is  also  a  contribution  to  from  the  current  in  the  upper  sheet.  The  computations 
for  this  integral  are  the  same  as  for  the  one  already  done  if  we  set  m  =  0,  except  that  the  sign 
is  opposite  because  the  currents  in  the  upper  stripline  are  the  negative  of  the  currents  in  the 
lower  stripline.  If  m  is  zero,  \  becomes  merely  k/b. 

For  the  upper  stripline  currents,  the  expression  corresponding  to  (B-2)  is 


oo 


1 

cb 


Z 


n=0 


A  H  (b) 
n  nq 


(B-3) 


H 


u 

nq 


(b) 


ps  (-i)^~q  gdKb 
4  K(Q^  -  N^) 


{2Q  Co(K)  sin  (3^  K) 


(b) 

qq 


4 


£ 


,  ,  NK ,  XT  •  /  QK ,  ,  g 

+  (Q  sin(— )  -  N  sin(^)  ]  |j 

g  bdK  (|j.  So(K)  I(K)  - 


K 


So(K)}  , 

Co(K)  sin(^^)  ] 
_ 2 _ 


Therefore,  we  find  that  after  dividing  out  the  sinusoidal  dependence  on  y,  multiplying  by 
cos(Qbx/g),  and  integrating  from  —  W/2  to  W/2,  Eq.  (7)  reduces  to 
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4-  y  A  H  (b) 
be  ^  n  nq  ’ 

n=0 

t,  /VI  ^  r®  (-1)""'’ gdK  .,Q  +  N„.  -  /Q-N^. 

H  (b)  =  4  \  2  2  '  ~r  °c 

<Jo  K\  (Q  -  N  )  S  ® 

+  [N  sin(^)  -Q  sin(^^)]  1^  Dg} 

H  ^(b)  =  4  r®  £  [b\^  So(K)  -  K  So(b\)]  I(K) 

^  Kx'"  1^^ 

+  [^  Co{h\)  -h\^  Co(K)j|  . 

We  now  consider  (V  X  J)  •  1^  (pc/s)  and  execute  the  same  operations  applied  to  (7) 
From  (9), 

OO 

•(Vxj)=b^y  A  (H  +  g  )  cos(— )  sin(by) 
sz  s^  n  g  N  e 


n=0 


Dividing  by  sin  (by),  multiplying  by  cos(Qbx/g),  and  integrating: 


bpc 

~s~ 


V  A  /N  ,  gv  /Nbx.  /Qbx.  ,  bpc  g  .  /Q  .  g\ 

^  n'g  ^  J_„/2  g  g  s  7  p'g  Q' 


g/2 


n=0 


-g/2 


Therefore,  with  H^^(b)  given  by  (B-4): 


Z  An  V) 

n=0 


(B-4) 


(B-5) 
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